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ABSTRACT 



Aims. We reassess the properties of the ICM at large radii in the galaxy cluster PKS 0745-191 in light of recent Suzaku 
measurements. 

Methods. We analyze an archival 10.5 ksec _R0S'j4r/PSPC observation to extract the surface-brightness profile of PKS 
0745-191 and infer the deprojected density profile. We then compare the ROSAT surface-brightness profile with the 
Suzaku result. We perform a mass analysis by combining the ROSAT density profile and the published temperature 
profiles from different instruments. 

Results. We find that the ROSAT surface-brightness profile is statistically inconsistent (7.7(t) with the Suzaku result 
around and beyond the value of r2oo estimated by Suzaku. We argue that, thanks to its large field of view and low 
background, _R0S'j4r/PSPC is to the present day the most sensitive instrument for studying low surface-brightness 
X-ray emission in the 0.4-2.0 keV band. We also note that the Suzaku temperature and mass profiles are at odds with 
the results from at least two other satellites {XMM-Newton and Swift). 

Conclusions. The difference in surface brightness between ROSAT and Suzaku is most likely explained by the existence 
of additional foreground components at the low Galactic latitude of the source, which were not taken into account 
in the Suzaku background modeling. In light of our mass analysis, we conclude that any estimate of the fraction of 
the virial radius reached by X-ray measurements is affected by systematic errors of the order of 25%. As a result, the 
properties of the ICM at the virial radius are still uncertain, and the Suzaku results should be considered with caution. 
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1. Introduction 

Observing the outskirts of galaxy clusters is important for 
understanding the formation processes of large-scale struc- 
tures. For instance, studying the intracluster gas around 
the virial radius gives an opportunity to measure the tran- 
sition between the gravitationally-bound gas of clusters and 
the infalling material from large-scale structures. Tracing 
the gas out to the virial radius is also important for cal- 
ibrating the X-ray ma ss measure ments, which are impor- 
tant to cosmology fe.g.. lVoitll2005[ ). However, measuring the 
properties of the intracluster medium in cluster outskirts 
is difhcult because of the low surface brightness of these 
regions. An instrument with a low internal background is 
required to measure the source emission at high radii. In 
the past few years, the Suzaku satellite achieved a break- 
through in this domain, performing measurements of the 
ther modynamical properties of the ICM out to the virial ra- 
dius (iBautz et al.ll2009HReiprich et "alll2009t E oshino et al.l 



l2010HKawaharada et al.ll2010[ ). 



PKS 0745-191 {z = 0.1028) is a ve ry luminous (Lx 
3x10^^^ ergs s~^ in the 2-10 keV band. lArnaud et al.|[T987t 
I Allen et al.lll996h . cool-core cluster located in the vicinity 
of the Galactic plane (b — +3°). From a S uzaku/XlS obser- 
vation of the cluster, iGeorge et al.l ()2009D (hereafter, G09) 
measured the cluster emission out to ~ 1.5r2oo, and deter- 
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mined a value of r2oo = 1-7 Mpc (15.2 arcmin) for the virial 
radius. Surprisingly, the authors noted a flattening of the 
density and entropy profiles around r2oo, at variance with 
resul t s from cosmological sim ulations fe.g.. lRoncarelli et al.l 
120061 : iTozzi Sz Normanll2001l ) . This and similar results from 
other authors inspired a large amount of t heoretical work 
(e.g., iLapi et al.l 120101: iNagai fc Laull201lD . invoking sev- 
eral mechanisms (e.g., non-thermal pressure support, gas 
clumping) to reconcile simulations and observations. An 
independent confirmation of this result would therefore be 
very important to our understanding of cluster outskirts. 

In this paper, we present the analysis of an archival 
ROSAT Position Sensitive Proportional Counter (PSPC) 
observation of PKS 0745-191, with the aim of confirming 
the result of G09. Although the PSPC could not measure 
temperatures because of its limited bandpass (0.1-2.4 keV), 
its low instrumental background and large field of view (~2 
square degrees) made it an excellent tool for the study of 
low surface-brightness reg ions such as the outer regions of 
galax y clusters (see e.g., IVikhlinin et al.l [l999l : iNeumaniil 
I2005D . We also perform a mass analysis using the PSPC 
density profile and the temperature measurements from 
various other X-ray satellites, and compare the results with 
the measurements of G09. The paper is organized as fol- 
lows. In Sect. 121 we describe the data analysis procedure. 
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Fig. 1. Background-subtracted (NXB+SSX), exposure- 
corrected, adaptively smoothed ROSAT /PSPC image of 
PKS 0745-191. The black areas indicate the position of 
masked point sources. The white circle shows the approxi- 
mate location of r2oo as estimated by G09. 
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Fig. 2. Fit of the SB profile in the 25-50 arcmin radial 
range by a constant to estimate the cosmic background 
component (blue). For comparison, the NXB+SSX profile 
is shown in red. The bottom panel shows the deviations of 
the data points to the model in units of a. 



We present our results for the density profile of the cluster 
in Sect. [31 and discuss them in Sect. 01 

Throughout the paper, we assume a ACDM cosmology 
with nm = 0.3, Qa = 0.7, and Hq = 70 km s^^ Mpc"!. 

2. Data analysis 

2.1. Reduction 

PKS 0745-191 was the target of a pointed ROSAT /PSPC 
observation on October 15, 1993 (observation ID 
RP800623N00) for a total of 10.5 ksec. We reduced the 
data usi ng the ROSAT Exten ded Source Analysis Software 
(ESAS, ISnowden et al.l |1994[ ). To eliminate flaring peri- 
ods, we extracted a light curve from the raw event files 
and rejected all time periods where the Master Veto count 
rate was above 220 counts s~^. We then used the ESAS 
task ao to create a model of the scattered solar X-ray 
background fSSX. [Snowden &: Frevber3ll993() . and gener- 
ated a p a,rticle background rnqdel using the cast-pa rt ex- 
ecutable (jSnowden et al.lll992l: [Plucinskv et al.lll993D . The 
total model for the non X-ray background (NXB) and the 
SSX was then inferred. 

A counts image in the 0.4-2.0 keV band was then ex- 
tracted from the cleaned event file and the correspond- 
ing exposure map was created using the task cast_exp. 
Point sources were then detected from the image (using 
the program detect) and a point source mask was gener- 
ated to excise the corresponding regions. As a result, a 
background-subtracted, exposure-corrected image was cre- 
ated and adaptively smoothed. In Fig. [I] we show the re- 
sulting image together with the position of r2oo estimated 
by G09. Point sources have been masked from the image. 

2.2. Density profile 

To measure the density profile of the cluster, we first ex- 
tracted the surface-brightness (SB) profile of the source fol- 



lowing the procedure described in lEckert et al.l (|201lt) . We 
used the total counts image and the exposure map to ex- 
tract the mean count rate in concentric annuli out to a 
radius of 50 arcmin with a bin size of 1 arcmin. The NXB 
and SSX profiles were also extracted using the same proce- 
dure from the background image and subtracted from the 
total profile. 

To subtract the cosmic background component, we fit- 
ted the SB profile in the 25-50 arcmin radial range (r > 
l-6?'2oo) with a constant (see Fig. [5]). The bins were grouped 
to ensure a minimum number of 200 counts per bin. The 
profile in this radial range is accurately-described by a con- 
stant (x^ = 18.6/21 d.o.f.). We verified that the back- 
ground determination does not change significantly when 
using a different radial range by repeating the same pro- 
cedure in the radial ranges 25-45, 30-45, 25-40, 30-50, and 
35-50 arcmin. In all cases, we find that the measured back- 
ground value is consistent to within la with the value ex- 
tracted from the total 25-50 arcmin range, which indicates 
that our background determination is stable. The resulting 
background value was then subtracted from the total pro- 
file, and the uncertainties in the background determination 
were propagated to the background-subtracted profile. In 
Fig. [21 we show the surface-brightness profile in the 25-50 
arcmin range together with the best-fit value for the sky 
background (blue) and the NXB-fSSX profile (red). We 
can see that at all radii the cosmic component is dominant 
relative to the NXB and the SSX. 

From the background-subtr acted profile, we used the 
procedure of iKriss et al.l (|1983[ ) to deproject the surface- 
brightness profile. We assumed the temperature profile of 
G09, and used XSPEC vl2.6.0 and the PSPC effective area 
to convert the PSPC count rate per unit volume into emis- 
sivity using the normalization of the MEKAL model 

^°™=4.[d^(n-z)]^y"^"^'^^' 
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Fig. 3. ROSAT /PSPC density profile. The error bars and 
upper limits were estimated from Monte Carlo simulations 
(see text). The dotted black line shows the approximate 
location of r2oo computed by G09. 



where ~ 1.2n^f. We note that in the 0.4-2.0 keV band 
the conversion factor is insensitive to temperature. For tem- 
peratures between 2 and 8 keV, the conversion from PSPC 
count rate to emissivity changes by at most 4%. Assuming 
spherical symmetry and a constant density within each 
shell, we deduced the density profile. 

2.3. Error estimate 

A crucial point for our analysis is the estimate of statis- 
tical uncertainties in the density profile. To compute the 
error bars in the density profile, we used a Monte Carlo ap- 
proach to analyzing the ROSAT image. Assuming Poisson 
statistics, we generated 10^ realizations of the PSPC counts 
profile, and performed the procedure described above to ob- 
tain 10* realizations of the density profile. The error bars 
and confidence intervals were then calculated from the dis- 
tribution of values in each density bin. When a negative 
background-subtracted profile was found, the density was 
set to 0. We then defined a 90% upper limit by the value 
for which 90% of the simulations give a result below this 
value. 



3. Results 

In Fig. H we show the resulting ROSAT /PSPC density 
profile. The dotted black line shows the estimate of r2oo 
from G09. Beyond r = 17', we do not detect any significant 
cluster emission, and set an upper limit to the density in 
the 17 - 25' range of ni7-25' < 4.2 x IQ-^ cm'^ (90% 
confidence level). 

To compare the ROSAT and Suza ku profiles, we used a 
metho d similar to the one reported in lBowver fc VikhlininI 
(l200i). More specifically, we converted the projected 
Suzaku profile into the MEKAL normalization (Eq. [ij, 
and folded it through the response of ROSAT /PSPC to 
compute the count rate that should be seen by PSPC. 
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Fig. 4. Background-subtracted (cosmic-f NXB) 

ROSAT /PSPC surface-brightness profile (red), com- 
pared to the Suzaku profile from G09 folded through the 
PSPC response (black). 

The resulting profile (black) , compared to the background- 
subtracted PSPC profile with the same binning (red), is 
shown in Fig.|4l While in the inner 4 annuli the profiles dif- 
fer by only ~ 11% (the difference in the innermost 2 bins 
being explained by the broader Suzaku PSF), a clear dis- 
crepancy between the ROSAT and Suzaku profiles is found 
beyond 13.5'. In the 13.5-18.5 arcmin bin, the Suzaku pro- 
file exceeds the ROSAT data point by a factor of three, 
while in the 18.5-24 arcmin bin our 90% upper limit lies a 
factor of 2.1 below the Suzaku detection. The profiles at the 
largest radii are discrepant at more than 7.7cr. This result 
is stable in terms of the background level. We indeed reach 
the same conclusion when we use the lowest allowed value 
for the background instead of the mean value. 

4. Discussion 

As shown above, our ROSAT /PSPC density profile is sta- 
tistically inconsistent with the Suzaku profile from G09 at 
high significance. We discuss here several possibilities to 
reconcile the two results. 

1. Cross- calibration: As found in the case of 
ROSAT/PSPC a nd XMM-Newtor i /M0S2 (see 
Appendix A of lEckert et all l201lt) . there could 
be a cross-calibration discrepancy between PSPC and 
Suzaku /I^IS. However, as shown in Fig. 2] below 13.5 
arcmin the projected profiles agree to within 11%, 
while at the largest radii the ROSAT data points are 
below the Suzaku measurements by a factor of 2.1 — 3.0. 
Thus, although a cross-calibration problem may ex- 
ist, it certainly cannot explain the observed discrepancy. 

2. Absorption: Because of its softer bandpass, ROSAT 
is more affected than Suzaku by absorption, which is 
large at the position of PKS 0745-191 [Nh = 4.2 x lO^^ 
cm~^). If the Galactic absorption is significantly larger 
in the outskirts of the cluster {Nh ^ 10^^ cm~^) than 
in the center, the flux observed in the Suzaku data can 
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be suppressed in the ROSAT band. Variations in the 
Nh on these scales (~ 200 arcmin'^) should produce 
significant variations (by a factor of ^ 2.2) in the CXB 
intensity in the ROSAT band. However, an analysis 
of the background beyond 25 arcmin from the cluster 
center did not show any significant variations in the 
background level. It therefore seems very unlikely that 
these variations are present only along the line of sight 
of the cluster. 

3. Spatial coverage: Because of the smaller field of view 
(FOV) of Suzaku, the G09 result was extracted in 
a rather small region (see Fig. 1 of G09). If the gas 
density v aried significantly at large radii (as in the case 
of A1795. lBautz et al.ll2009l ). the regions selected in the 
Suzaku observation might have a higher density, hence 
would not be representative of the mean density around 
the virial radius. However, extracting the surface- 
brightness profile in exactly the same region as G09, we 
do not find any significant difference from the mean pro- 
file. In this specific region, the discrepancy with Suzaku 
is still statistically significant {5.2a in the last two bins). 

4. Background determination: Because of the very low sur- 
face brightness of the cluster emission around r2oo , the 
background determination is important when measur- 
ing the physical properties of the gas in these regions. 
Thanks to its much wider FOV, ROSAT has a clear ad- 
vantage with respect to Suzaku. The large ROSAT FOV 
allows us to perform a local background measurement, 
while G09 had to rely on background modeling. At the 
low Galactic latitudes where PKS 0745-191 is located, 
the foreground emission associated with our galaxy is 
known to be stronger and more variable than at higher 
latitudes, making background characterization a more 
complicated business. Thanks to its higher spatial res- 
olution and larger FOV, ROSAT is able to more accu- 
rately measure the total sky background emission, and 
resolve a larger fraction of sources (see Fig. [T]) , so the 
systematic uncertainties affecting the ROSAT measure- 
ments are smaller than those of Suzaku. 

An additional background component is probably caus- 
ing the discrepancy. G09 used an earlier observation of 
the Lockman hole to estimate the cosmic background 
component, arguing that, since the total cosmic back- 
ground level measured by ROSAT in the 1-2 keV band 
at the position of PKS 0745-191 is similar to the one for 
the Lockman hole, the same model could be used for the 
sky background in the two observations. However, since 
the extragalactic background at the low Galactic lati- 
tude of PKS 0745-191 is lower because of the larger ab- 
sorption, another component, most likely Galactic, must 
be responsible for com pensating t h e dec rease in the ex- 
tragalactic emission. Masui et al.l (|2009[ ) presented the 
Suzaku spectrum of an empty field located only 3° away 
from PKS 0745-191. The authors found clear evidence 
of another sky component at a temperature of '~ 0.8 
keV, with excess emission in the 1-1.2 keV range sug- 
gesting that even higher temperature emission exists. 
This implies that the spectrum is qualitatively unlike 
empty-field spectra extracted at high Galactic latitude. 
Since G09 neglected these components, their results are 
likely to be affected by an incorrect background model, 
implying that clusters located at low Galactic latitude 



(or in any region with complicated soft emission, such 
as the North polar spur) are not ideal targets for the 
study of cluster outskirts. 

In low-SB regions, the key feature for an instru- 
ment is the ratio of b ackground level to effective area 
(|Ettori fc Molendi|[20T0l ). In this respect, ROSAT /PSPC 
is a more sensitive instrument than SuzakufKIS. From the 
total background level (see Fig. [2]) we get a value of 4 x 10~^ 
counts s~^ keV~^ arcmin"^ cm~^ for this quantity (0.4-2.0 
keV band), which is almos t a factor of two low er than for 
Suzaku /XIS (see Fig. 5 of 'Mitsud a et aI1l2007^ . 

As stated in G09, the mass parameters derived from the 
Suzaku data appear to be inconsistent with other results. 
We performed a mass analysis using the ROSAT density 
profile and the temperature profiles from various satellites 
{Suzaku, XMM, BeppoSAX and Swift, see Appendix A). 
Interestingly, we see that all satellites except Suzaku find 
a value of r2oo that is larger than 2 Mpc (~ 25% larger 
than the value estimated by G09), corresponding to a dif- 
ference of a factor of 2 in the virial mass. For XMM and 
Swift measurements the difference in r2oo is statistically 
highly significant, at more than 6 and 5(T, respectively (see 
Appendix A for details). 

Since there is no compelling reason to prefer one set 
of measurements with respect to another, we are forced 
to conclude that the scale radius in PKS 0745-191 is cur- 
rently affected by a systematic indetermination of roughly 
25%. It goes almost without saying that any estimate of 
the fraction of the virial radius reached by X-ray mea- 
surements will have to take this indetermination into ac- 
count. This is no small issue. Depending on whether we 
assume the Suzaku or the XMM / Swift virial radius we 
have that the PSPC surface brightness measurements ex- 
tend to 1.03r2oo or 0.8r2oo- Finally, although we cannot 
prefer one estimate of r2oo to another, there are argu- 
ments favoring the XMM / Swift measurement over the 
Suzaku one. First of all, we have two independent mea- 
surements that are consistent with one another and in- 
consistent with the Suzaku one. Secondly, the Swift and 
XMM profiles are consistent with me an cluster temperature 
profil es measur ed with BeppoSAX (iDe Grandi fc Molendil 
'200 21). XMM (iLeccardi fc Molendil I2008D. and Chandra 



(Vi khlinin et al 



2006f). and with predictions from simula- 



tions (e.g.. iRoncarelli et ani2006|l . This is not the case for 
the Suzaku profile, which shows a much more rapid and 
significant decline. This is an important point since, as dis- 
cussed in the Appendix, the difference in r2oo follows from 
the difference in the shape of the radial temperature pro- 
files. 



5. Conclusion 

We have reported our analysis of an archival 
ROSAT /PSPC observation of the galaxy cluster PKS 
0745-191. We have found that the surface-brightness profile 
extracted from PSPC data is statistically inconsistent with 
the Suzaku result from G09 at high significance {7.7a). 
At large radii (> 13.5 arcmin), the predicted count rate 
exceeds the PSPC data by a factor of 2.1 — 3. This 
difference is most likely due to a problem in the modeling 
of the background in the Suzaku measurement, which is 
difficult for a narrow-field instrument at the low Galactic 
latitude of the source. Thanks to its larger FOV and 
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higher spatial resolution, ROSAT is able to measure the 
background more accurately than Suzaku. 

We have also shown that the rapidly declining temper- 
ature profile measured by Suzaku leads to total mass and 
7'2oo estimates at variance with those derived from Swift 
and XMM temperature profiles. In the absence of com- 
pelling proof favoring one measurement over another, we 
conservatively conclude that any estimate of the fraction 
of the virial radius reached by X-ray measures is affected 
by a systematic error of about 25% associated with the 
indetermination in r2oo- In conclusion, we stress that the 
current observational knowledge of the ICM properties of 
PKS 0745-191 at large radii is still uncertain, thus the re- 
sults presented by G09 should be considered with care. 
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Appendix A: Mass analysis with different datasets 

We perforin e d a m ass analysis adopting method 2 of 
lEttori et al.l (|2010t) (and references therein) using the 
ROSAT surface brightness profile (this work). For the 
temperature profile, we adopted the published tempera- 
tu re profiles of four different instr uments: XMM-Newton 
(S nowden et alj [20081. B eppoSAX (iDe Grandi fc Molendil 
200'3), Swift /XRT ( Moretti et al.ll2010f ). and the Chandra 
and Suzaku temperature profile of G09 (their Fig. 5, i.e. ex- 
cluding the Chandra data points within 70 kpc and the 
innermost and outermost Suzaku data points) . In Fig. lA.ll 
we plot the temperature profiles for PKS 0745-191 within 
the inner 15 arcmin from the center of the cluster, obtained 
by various satellites with very different instruments (gas 
proportional counters or CCDs) in very different orbits (low 
Earth orbits (LEO) with very low particle background such 
as BeppoSAX and Swift /XRT, or highly elliptical orbits 
with high background such as XMM). There is a clear dis- 
crepancy between the behavior of the Suzaku profile com- 
pared to all the other satellites. 

In Table lA.li we give the best-fit parameters for the 
NFW profile estimated using the different datasets. We can 
see that XMM-Newton, BeppoSAX, and Swift data give 
results in good agreement, and point towards a value of 
r2oo larger than 2 Mpc. In contrast, if we use the Chandra 
and Suzaku profile of G09 we obtain results in agree- 
ment with theirs, r2oo = 1-88 ± 0.03 Mpc and Af2oo = 
(8.32±0.41) X 10"Mo, notwithstanding the different meth- 
ods adopted in the mass analysis, for example the use of 
the ROSAT density profile and the simplifying assump- 
tions of constant density within each shell and estimate of 
the volume defined by G09. G09 correctly points out the 
problems of the NFW fits applied to the Chandra data be- 
cause the data do not extend to the NFW scale radius, r^. 
The mass fits reported in Table ED all use the ROSAT 
surface brightness profile extending to 1927 kpc and the 
spectroscopic temperature data out to the radius listed as 
r^sp for the various satellites. A necessary ingredient for a 
reliable measurement of the NFW parameters is that the 
fitted value of rs lies well within the outer radius of the 
X-ray data, which is fulfilled by all the adopted datasets. 
Between XMM-Newton and Suzaku, the discrepancy is at 
the 6a and 4.8(t level for r2oo and M200, respectively. The 
results are at odds at the level of 5.6a and 4.6cr between 
Swift and Suzaku, and 1.5ct and 1.2cr between BeppoSAX 
and Suzaku. 

These results clearly highlight the importance of the 
Suzaku temperature profile for the rather low values of r2oo 
and M200 obtained in G09 compared to those obtained from 
XMM-Newton, BeppoSAX, and Swift data sets. The dis- 
crepancy can be quantified by comparing the value for the 
slope of the power-law fit to the temperature profile exclud- 
ing the core as reported by G09, —0.94 ± 0.06, to the mean 
values of observed temperature profiles beyond 0.2ri8o as 
reported by Leccardi & Molendi (2008), -0.31 ± 0.02, i.e. 
a 10a deviation from the mean. The inconsistency with the 
milder decline found in hydrodynamical simulations (e.g., 
iRoncarelli eralll2006[ ) was also noted by G09. 



6 



Eckert, D. et al.: X-ray observations of PKS 0745-191 at the virial radius: Are we there yet? 



Table A.l. Best-fit parameters of the NFW profile to the ROSAT density and several published temperature profiles (see 
Fig. I A. l]) . The parameters are the concentration (c2oo)j NEW scale radius (rj) in kpc, r2oo in kpc, and gravitational mass 
(Af2oo) in units of 10^'* M0. The last parameter, Vxsp, denotes the maximum radius at which temperature measurements 
could be performed, in kpc. 





ROSAT+XMM 


ROSAT+Suzaku 


ROSAT+SAX 


ROSAT+Swift 


C200 


3.45 ±0.25 


9.05 ±0.16 


4.89 ± 0.89 


3.52 ±0.26 




697 ± 77 


207 ±2 


434 ± 137 


689 ± 76 


^200 


2405 ± 82 


1878 ± 33 


2124 ± 160 


2429 ± 92 


M20O 


17.48 ± 1.85 


8.32 ±0.41 


12.03 ± 3.04 


18.00 ± 2.05 


xVd.o.f. 


57/23 


42/20 


27/18 


61/19 


fxsp 


1587 


2097 


1360 


1247 



-4 



XMM 
SUZAKU 
SAX 
SWIFT 



10 



15 



r (arcmin) 



Fig. A.l. Temperature profiles of PKS 0745-191 mea- 
sured by Su zaku (red tr iangles, G09), XMM-Newton (black 
circles, .Snowden et all [200S) . Swift /XRT (blue squares, 



Moretti et all I2OIOD and BevvoSAX ( green diamonds, 



De Grandi fc Molendil[2002h 



